Extremal entanglement and mixedness in continuous variable systems 
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We investigate the relationship between mixedness and entanglement for Gaussian states of continuous vari- 
able systems. We introduce generalized entropies based on Schatten p-norms to quantify the mixedness of a 
state, and derive their explicit expressions in terms of symplectic spectra. We compare the hierarchies of mixed- 
ness provided by such measures with the one provided by the purity (defined as tr q 2 for the state q) for generic 
n-mode states. We then review the analysis proving the existence of both maximally and minimally entangled 
states at given global and marginal purities, with the entanglement quantified by the logarithmic negativity. 
Based on these results, we extend such an analysis to generalized entropies, introducing and fully characterizing 
maximally and minimally entangled states for given global and local generalized entropies. We compare the 
different roles played by the purity and by the generalized p-entropies in quantifying the entanglement and the 
mixedness of continuous variable systems. We introduce the concept of average logarithmic negativity, showing 
that it allows a reliable quantitative estimate of continuous variable entanglement by direct measurements of 
global and marginal generalized p-entropies. 
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I. INTRODUCTION 

The degree of entanglement (i.e. the contents in quantum 
correlations) as well as the degree of mixedness (i.e. "the 
amount" by which a state fails to be pure) are among the cru- 
cial features of quantum states from the point of view of quan- 
tum information theory. Indeed, the search for proper analyt- 
ical ways to quantify such features for general (mixed) quan- 
tum states cannot be yet considered accomplished. In view of 
such considerations, it is clear that the full understanding of 
the relationships between the quantum correlations contained 
in a bipartite state and the global and local (i.e. referring to the 
reduced states of the two subsystems) degrees of mixedness of 
the state, would be desirable. In particular, it would be a rel- 
evant step towards the clarification of the nature of quantum 
correlations and, possibly, of the distinction between quan- 
tum and classical correlations of mixed states, which remains 
an open issue |l||. A simple question one can raise in such 
a context is the investigation of the properties of extremally 
entangled states for a given degree of mixedness. 

Let us mention that, as for two-qubit systems, the notion 
of maximally entangled states at fixed mixedness (MEMS) 
was originally introduced by Ishizaka and Hiroshima ]2|]. The 
discovery of such states spurred several theoretical works 
01, aimed at exploring the relations between different mea- 
sures of entanglement and mixedness |4|| (strictly related to 
the question of the ordering of these different measures j5|). 
Moreover, maximally entangled states for given local (or 
"marginal") mixednesses have been recently introduced and 
analyzed in detail in the context of qubit systems @|. On 
the experimental side, much attention has been devoted to 
exploring the two-qubit Hilbert space in optical settings |7], 
while the experimental realization of MEMS has been re- 
cently demonstrated 1 8 ] . 

Because of the great current interest in continuous variable 
(CV) quantum information H G3 El 0, 

the extension of 

such analyses to infinite dimensional systems is higly desir- 



able. In the present work, we introduce and study in detail 
extremally entangled mixed Gaussian states of infinite dimen- 
sional Hilbert spaces for fixed global and marginal generalized 
entropies, significantly generalizing the results derived earlier 
in Ref. 1 13], where the existence of maximally and minimally 
entangled mixed Gaussian states at given global and marginal 
purities was first discovered. In the present paper we will 
make use of a hierarchy of generalized entropies, based on 
the Schatten p-norms, to quantify mixedness and character- 
ize extremal entanglement in continuous varibale (CV) sys- 
tems and investigate several related subjects, like the ordering 
of such different entropic measures and the relations between 
EPR (Einstein-Podolsky-Rosen) correlations and symplectic 
spectra. The crucial starting point of our analysis is the obser- 
vation that the existence of infinitely entangled states fl4HT^l 
prevents maximally entangled Gaussian states from being de- 
fined as states with maximal, finite, logarithmic negativity, 
even at fixed global mixedness 116|. However, we will show 
that fixing the global and the local mixednesses allows to de- 
fine unambiguously both maximally and minimally entangled 
mixed Gaussian states. This relevant and somehow surprising 
result - there is no analog of Gaussian minimally entangled 
states in finite dimensional systems - turns out to have an ex- 
perimental interest as well 11131 fl7ll . 

The paper is structured as follows. In Sec.||l]we briefly re- 
view the basic notation and the general properties of Gaussian 
states. In Sec. [Hi] we introduce the hierarchy of generalized 
Schatten p-norms and entropies, and we extensively discuss 
the problem of the ordering of different entropic measures 
for states with an arbitrary number of modes. In Sec. UVI we 
review the state of the art on the existing, computable mea- 
sures of entanglement for two-mode Gaussian states, while 
in Sec. [V] we present a heuristic argument relating EPR cor- 
relations to symplectic spectra. In Sec. I VII we introduce a 
parametrization of two-mode Gaussian states in terms of sym- 
plectic invariants endowed with a direct physical interpreta- 
tion. Exploiting these results, in Sec. lVUI we define maximally 
and minimally entangled states for given global and marginal 
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purities and present some experimental situations in which 
these states occur. In Sec. lVllll we generalize the concept of 
extremal entanglement in CV systems by introducing max- 
imally and minimally entangled states for given global and 
local generalized entropies of arbitrary order, and we present 
an extensive study of their properties. We find out, somehow 
surprisingly, that maximally and minimally entangled states 
can interchange their roles for certain ranges of values of the 
global and marginal p-entropies; moreover, we observe that 
with increasing p the generalized entropies, while carrying in 
general less information on a quantum state, provide a more 
accurate quantification of its entanglement. In Sec. lIXI we in- 
troduce and study the concept of "average logarithmic nega- 
tivity", showing that this quantity provides an excellent quan- 
titative estimate of CV entanglement based only on the knowl- 
edge of global and local entropies. Finally, in SecQ^we sum- 
marize our results and discuss future perspectives. 



n. GAUSSIAN STATES: GENERAL OVERVIEW 

Let us consider a CV system, described by an Hilbert space 
H = <S)"=i Hi resulting from the tensor product of the infi- 
nite dimensional Fock spaces Hi's. Let a; be the annihilation 
operator acting on Hi, andi; = (a.i+a\) andp; = (a,i — a\)/i 
be the related quadrature phase operators. The corresponding 
phase space variables will be denoted by Xi and p,-. Let us 
group together the operators ii and Pi in a vector of opera- 
tors X — (xi,pi, . . . ,x n ,p n ). The canonical commutation 
relations for the X;'s are encoded in the symplectic form fi 



with f2 



[Xi,Xj] = 2iQi 



0< 

i=l 



1 

-1 



The set of Gaussian states is, by definition, the set of states 
with Gaussian characteristic functions and quasi-probability 
distributions. Therefore a Gaussian state g is completely 
characterized by its first and second statistical moments, 
which form, respectively, the vector of first moments X = 

((X±), (Xi), . . . , (X n ), (X n )] and the covariance matrix of 

elements a 



1 



(1) 



where, for any observable 6, the expectation value (6) = 
Tr (go). First statistical moments can be arbitrarily adjusted 
by local unitary operations, which cannot affect any property 
related to entanglement or mixedness. Therefore they will be 
unimportant to our aims and we will set them to in the fol- 
lowing, without any loss of generality. Throughout the paper, 
<7 will stand for the covariance matrix of the Gaussian state g. 

Let us now consider the hermitian operator y = YX T , 
where Y E R 2 ™ is an arbitrary real 2rt-dimensional row vec- 
tor. Positivity of g imposes Tr(gy 2 ) > 0, which can be sim- 
ply recast in terms of second moments as YtY t > 0, with 



Tij = (XiXj). From this relation, exploiting the canonical 
commutation relations and recalling definition Q and the ar- 
bitrarity of Y, the Heisenberg uncertainty principle can be re- 
cast in the form 



(2) 



Inequality (|2j is the necessary and sufficient constraint a has 
to fulfill to be a bona fide covariance matrix 

fiHU. We 

mention that such a constraint implies cr > 0. 

In the following, we will make use of the Wigner quasi- 
probability representation W(xi,pi) defined, for any density 
matrix, as the Fourier transform of the symmetrically ordered 
characteristic function l20ll . In Wigner phase space picture, 
the tensor product H = &)Hi of the Hilbert spaces Hi's of 
the n modes results in the direct sum r — (J) of the phase 
spaces IYs. In general, as a consequence of the Stone-von 
Neumann theorem, any symplectic transformation acting on 
the phase space T corresponds to a unitary operator acting on 
the Hilbert space H, through the so called metaplectic rep- 
resentation. Such a unitary operator is generated by terms 
of the second order in the field operators fl9il . In what fol- 
lows we will refer to a transformation Si = (J) Si, with each 
Si G Sp(2,j&) acting on Ti, as to a "local symplectic opera- 
tion". The corresponding unitary transformation is the "local 
unitary transformation" Ui = Ui, with each Ui acting on 
Hi. 

The Wigner function of a Gaussian state can be written as 
follows in terms of phase space quadrature variables 



-hXc 



W(X) 



7rvDet (7 



(3) 



where X stands for the vector (x\,p\, . . . , x n ,p n ) £ T. 

Finally let us recall that, due to Williamson theorem I2TI1 . 
the covariance matrix of a n-mode Gaussian state can always 
be written as fl8il 



a = S T uS , 



(4) 



where S £ Sp^n,R) ar, d v is the covariance matrix 

v= diag(^i,^i, ...,u n ,v n ), (5) 

corresponding to a tensor product of thermal states with diag- 
onal density matrix given by 



fc=0 



1 



v, 



|*>«<*| 



(6) 



\k)i being the k-th number state of the Fock space Hi. The 
dual (Hilbert space) formulation of Eq. 10} then reads: g = 
U* g® U, for some unitary U. 

The quantities Ui's form the symplectic spectrum of the co- 
variance matrix a and can be computed as the eigenvalues of 
the matrix |zf2er|. Such eigenvalues are in fact invariant under 
the action of symplectic transformations on the matrix cr. 
The symplectic eigenvalues vi encode essential informations 
on the Gaussian state cr and provide powerful, simple ways to 
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express its fundamental properties. For instance, let us con- 
sider the Heisenberg uncertainty relation 0. Since S is sym- 
plectic, one has S^^CtS" 1 = fl, so that inequality (|2j is 
equivalent to v + iSl > 0. In terms of the symplectic eigen- 
values Vi the uncertainty relation then simply reads 



Vi > 1 



(7) 



We can, without loss of generality, rearrange the modes of a 
?i-mode state such that the corresponding symplectic eigen- 
values are sorted in ascending order 

V\ < V 2 < ■ ■ ■ < Vn-l < V-a ■ 

With this notation, the uncertainty relation reduces to v\ > 1. 
We remark that the full saturation of Heisenberg uncertainty 
principle can only be achieved by pure n-mode Gaussian 
states, for which vi = 1 Vi = 1, . . . , n. Instead, mixed 
states such that = 1 and v^t- > 1, with 1 < k < n, 
only partially saturate the uncertainty principle, with partial 
saturation becoming weaker with decreasing k. Such states 
are minimum uncertainty mixed Gaussian states in the sense 
that the phase quadrature operators of the first k modes satisfy 
the Heisenberg minimal uncertainty, while for the remaining 
n — k modes the state indeed contains some additional ther- 
mal and/or Schrodinger-like correlations which are responsi- 
ble for the global mixedness of the state. 



A. Two-mode states 

In the present work we will mainly deal with two-mode 
Gaussian states. Here, we will thus briefly review their rele- 
vant properties and specify some further notations. 

The expression of the two-mode covariance matrix er in 
terms of the three 2x2 matrices a, (3, 7 will be useful 



a 7 

7 T (3 



(8) 



It is well known that for any two-mode covariance matrix er 
there exists a local symplectic operation Si = Sj® S2 which 
takes er to the so called standard form cr s f 122.12311 



STer Si 



a c+ 

a c_ 

c+ b 

c_ b 



(9) 



States whose standard form fulfills a = b are said to be sym- 
metric. Let us recall that any pure state is symmetric and 
fulfills c + = — c_ = y/a 2 — 1. The correlations a, b, c+, 
and c_ are determined by the four local symplectic invariants 
Deter = (ab - c\){ab - <?_), Deta = a 2 , Det/3 = b 2 , 
Det7 = c + c_. Therefore, the standard form corresponding 
to any covariance matrix is unique. 

The uncertainty principle Ineq. (0 can be recast as a con- 
straint on the S'p(4 i r) invariants Deter and A(er) = Deta + 
Det/3 + 2Det7: 



A(cr) < 1 + Deter 



Let us mention that, as it is evident from Eq. me condition 
er > implies 



ab - 



4 > 



(11) 



The symplectic eigenvalues of a two-mode Gaussian state 
will be named v_ and v + , with v_ < v + in general. With 
such an ordering, the Heisenberg uncertainty relation Eq. @ 
becomes 



> 1 



(12) 



Full saturation v- = v+ = 1 yields the standard pure Gaus- 
sian states of Heisenberg minimum uncertainty; partial satura- 
tion V- = 1, > 1 defines the minimum-uncertainty mixed 
Gaussian states. A simple expression for the can be found 
in terms of the two Spu t «\ invariants j24ll25ll 



2*4 = A(er) T VA(^F 



4 Det er . 



In turn, Eq. Jl 3i yields immediately 
A(er) = v 2 _ + 



(13) 



(14) 



A subclass of Gaussian states that will play a relevant role 
in the following is constituted by the nonsymmetric two-mode 
squeezed thermal states. Let S r = exp(iraiQ,2 — irajaj) 
be the two mode squeezing operator with real squeezing pa- 
rameter r, and let £>® be a tensor product of thermal states 
with covariance matrix — 1.2 V - ® ^-2 V +, where v T is, 
as usual, the symplectic spectrum of the state. Then, a non- 
symmetric two-mode squeezed thermal state £„ 4|t . is defined 
as £„ it r = SrQuiS^,, corresponding to a standard form with 



a 
b 



V- cosh 



sinh 



c± = ±- 



+ v + sinh r 
+ v + cosh 2 r 

■ sinh 2r . 



(15) 



We mention that the peculiar entanglement properties of 
squeezed nonsymmetric thermal states have been recently an- 
alyzed 1 26]. In the symmetric instance (with v_ = v + = v) 
these states reduce to two-mode squeezed thermal states. The 
covariance matrices of these states are symmetric standard 
forms with 



a = v cosh 2r , c± = ±v sinh 2r . 



(16) 



(10) 



In the pure case, for which v = 1, one recovers the two- 
mode squeezed vacua. Notice that such states encompass all 
the standard forms associated to pure states: any two-mode 
Gaussian state can reduced to a squeezed vacuum by means 
of unitary local operations. 

Two-mode squeezed states (both thermal and pure) are en- 
dowed with remarkable properties related to entanglement 
ll27ll28ll29ll . Even the ideal, perfectly correlated original EPR 
state 1 30] can be seen indeed as a two-mode squeezed vacuum 
in the limit of infinite squeezing parameter r. The dynamical 
properties of the entanglement and the characterization of the 
decoherence of such states have been addressed in detail in 
several works ll3lll . We will also show, as a byproduct of the 
present work, the peculiar role they play as maximally entan- 
gled Gaussian states. 
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HI. MEASURES OF MIXEDNESS, DEGREE OF 
COHERENCE, AND ENTROPIC MEASURES 

The degree of mixedness of a quantum state g can be char- 
acterized completely by the knowledge of all the associated 
Schatten p-norms 13211 

\\g\\ p = {Tt\q\ p )p = (Tr qP)? , withp > 1. (17) 

In particular, the case p = 2 is directly related to the purity 
/i = Tr g 2 = (||f?||2) 2 13311 . The p-norms are multiplicative 
on tensor product states and thus determine a family of "gen- 
eralized entropies" S p 11341 13511 . defined as 



1 - Tr g p 

S p = — , p > 1. 

p- 1 



(18) 



These quantities have been introduced independently by M. 
J. Bastiaans in the context of quantum optics 1 34-], and by 
C. Tsallis in the context of statistical mechanics [35]. In the 
quantum arena, they can be interpreted both as quantifiers of 
the degree of mixedness of a state g by the amount of in- 
formation it lacks, and as measures of the overall degree of 
coherence of the state (the latter meaning was elucidated by 
Bastiaans in his analysis of the properties of partially coher- 
ent light). The quantity S% = 1 — /i = Sl, conjugate to 
the purity /i, is usually referred to as the linear entropy: it 
is a particularly important measure of mixedness, essentially 
because of the simplicity of its analytical expressions, which 
will become soon manifest. Finally, another important class 
of entropic measures includes the Renyi entropies 1 36] 



In Tr gP 
~1 



P 



p > 1. 



It can be easily shown that 



lim S„ = lim S p = — Tr (g\n g) = S\ 



(19) 



(20) 



so that also the Shannon-von Neumann entropy Sv can be de- 
fined in terms of p-norms. The quantity Sv is additive on ten- 
sor product states and provides a further convenient measure 
of mixedness of the quantum state g. 

It is easily seen that the generalized entropies S p 's range 
from for pure states to l/(p — 1) for completely mixed 
states with fully degenerate eigenspectra. Notice that Sv is 
infinite on infinitely mixed states, while Sl is normalized to 
1. We also mention that, in the asymptotic limit of arbitrary 
large p, the function Tr g p becomes a function only of the 
largest eigenvalue of g: more and more information about the 
state is discarded in such an estimate for the degree of purity; 
considering S p in the limit p — > oo yields a trivial constant 
null function, with no information at all about the state under 
exam. We also note that, for any given quantum state, S p is a 
monotonically decreasing function of p. 

Because of their unitarily invariant nature, the generalized 
purities Tr g p of generic n-mode Gaussian states can be sim- 
ply computed in terms of the symplectic eigenvalues Vi of a. 
In fact, a symplectic transformation acting on <x is embodied 



by a unitary (trace preserving) operator acting on g, so that 
Tr g p can be easily computed on the diagonal state v. One 
obtains 



TrgP = Y[g p (vi) 



(21) 



where 



2" 



[x + \)p -{x- iy 

A first consequence of Eq. Mil is that 

1 1 



Kq) = 



U. u i VDeter 



(22) 



Regardless of the number of modes, the purity of a Gaussian 
state is fully determined by the symplectic invariant Dot er 
alone. A second consequence of Eq. i21i is that, together with 
Eqs. Jl 8i and (I20> . it allows for the computation of the von 
Neumann entropy Sv of a Gaussian state g, yielding 



Sv(q) 



i=l 



(23) 



where 



x + 1 , ( x + 1 
In 



x — 1 f x — 1 
In 



Such an expression for the von Neumann entropy of a Gaus- 
sian state was first explicitly given in Ref. |H. Let us remark 
that, clearly, the symplectic spectrum of single mode Gaus- 
sian states, which consists of only one eigenvalue V\, is fully 
determined by the invariant Det cr — v\ Therefore, all the en- 
tropies Sp's (and Sv as well) are just increasing functions of 
Det a (i.e. of Sl) and induce the same hierarchy of mixedness 
on the set of one-mode Gaussian states. This is no longer true 
for multi-mode states, even for the relevant, simple instance 
of two-mode states. 

Here we aim to find extremal values of S p (for p ^ 2) for 
fixed Sl in the general n-mode Gaussian instance, in order 
to quantitatively compare the characterization of mixedness 
given by the different entropic measures. For simplicity, in 
calculations we will replace Sl with \i. In view of Eqs. (12 1 1 
and \22\ . the possible values taken by S p for a given /i are 
determined by 

(P-1)S P = 1 - ( lIgM ) g p ( 1 _ 1 ) , (24) 



with 1 < Si < 



(25) 



The last constraint on the n — 1 real auxiliary parameters Sj is 
a consequence of the uncertainty relation 0. We first focus 
on the instance p < 2, in which the function S p is concave 
with respect to any Si, for any value of the s^'s. Therefore 
its minimum with respect to, say, s„_i occurs at the bound- 
aries of the domain, for s„_i saturating inequality J25b . Since 
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S p takes the same value at the two extrema and exploiting 
g p (l) — 1, one has 



(p-1) min S p = 1-1 Y\ 9p( s 



9p 



1 



n7l— I 
i=i 



(26) 



Iterating this procedure for all the Si's leads eventually to the 
minimum value S v 
ply reads 



jpminiv) of <S P at given purity p,, which sim 



p-1 



p < 2. 



(27) 



For p < 2, the mixedness of the states with minimal general- 
ized entropies at given purity is therefore concentrated in one 
quadrature: the symplectic spectrum of such states is partially 
degenerate, with v\ = . . . = v n -\ = 1 and v n = I /p. 

The maximum value S p max (p) is achieved by states satis- 
fying the coupled trascendental equations 

1 \ , 1 



9p 



g'p(sj) 



ns 3 n s, 



9p( s ])g'p 



where all the products JT run over the index i from 1 to n 
and 



(28) 
1, 



s40) = 



- P 2P [(x + l) p_1 - (x- l)?" 1 ] 
[(x + 1)p - (x - 1)p} 2 



(29) 



It is promptly verified that the above two conditions are ful- 
filled by states with a completely degenerate symplectic spec- 
trum: v\ = . . . = v n = p~ x l n , yielding 



5. 



p max 



(a0 = 



1 - 9p (/« " ) 
~ 1 



p < 2. 



(30) 



The analysis that we carried out for p < 2 can be straight- 
forwardly extended to the limit p — > 1, yielding the extremal 
values of the von Neumann entropy for given purity /i of n- 
mode Gaussian states. Also in this case the states with max- 
imal Sv are those with a completely degenerate symplectic 
spectrum, while the states with minimal Sv are those with all 
the mixedness concentrated in one quadrature. The extremal 
values SvminilA and S V max{p) read 



Svmin(n) = f ( — 



(31) 



(32) 



The behaviors of the von Neumann and of the linear entropies 
for two-mode states are compared in Fig.^ 

The instance p > 2 can be treated in the same way, with 
the major difference that the function S p of Eq. (I24t is convex 
with respect to any Si for any value of the Si's. As a conse- 
quence we have an inversion of the previous expressions: for 
p > 2, the states with minimal 5 , pm i ra (/x) at given purity /x are 
those with a fully distributed symplectic spectrum, with 



1 



9p 



p-1 



p > 2. 



(33) 



5 
4 



0.2 



OA 



0.6 



0.8 



FIG. 1 : Plot of the curves (solid lines) of maximal and minimal von 
Neumann entropy at given linear entropy for two-mode Gaussian 
states. The density of states in the plane of entropies is represented 
as well, by plotting the distribution of 20000 randomly generated 
states (dots). 




FIG. 2: Plot of the absolute difference between the maximal and the 
minimal values of the generalized entropies S p at fixed linear en- 
tropy Sl for two-mode Gaussian states, as a function of p. Different 
curves correspond to different values of the linear entropy: Sl ~ 0.8 
(dashed line), Sl = 0.5 (continuous line), and Sl = 0.1 (dash- 
dotted line). 



On the other hand, the states with maximal S p ma2 
purity fi are those with a spectrum of the kind v\ 
v n -i = 1 and v n = l//i. Therefore 



-0) 



P-1 



p > 2. 



at given 



(34) 



Sp min | decreases 



As shown in Fig. |2 the distance \S pjn 
with increasing p. This is due to the fact that the quantity S p 
carries less information with increasing p, and the knowledge 
of p provides a more precise bound on the value of S p . 



IV. MEASURES OF ENTANGLEMENT 

We now review the main aspects of entanglement for CV 
systems and discuss some of the possible quantifications of 
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quantum correlations for Gaussian states. 

The necessary and sufficient separability criterion for a 
two-mode Gaussian state g has been shown to be positivity 
of the partially transposed state g (PPT criterion) 1 22] . In gen- 
eral, the partial transposition g of a bipartite quantum state g 
is defined as the result of the transposition performed on only 
one of the two subsystems in some given basis. Even though 
the resulting g does depend on the choice of the transposed 
subsystem and on the transposition basis, the statement g > 
is geometric, and is invariant under such choices 1 38 ] . It can 
be easily seen from the definition of W(X) that the action 
of partial transposition amounts, in phase space, to a mirror 
reflection of one of the four canonical variables. In terms of 
<Sp(2,R)ffiSp(2,R) invariants, this reduces to a sign flip in Det 7 



Det7 



-Dct 7 . 



Therefore the invariant A(er) is changed into A(er) = 
A(er) = Det a + Det (3 — 2 Det 7. Now, the symplectic 
eigenvalues v T of e> read 



\ 



A(er) =p JA(cr) 2 - 4Deter 



(35) 



The PPT criterion thus reduces to a simple inequality that 
must be satisfied by the smallest symplectic eigenvalue V- of 
the partially transposed state 



v-> 1, 



which is equivalent to 



A(er) < Deter + 1 



(36) 



(37) 



The above inequalities imply Det 7 = c + c_ < as a nec- 
essary condition for a two-mode Gaussian state to be entan- 
gled. The quantity encodes all the qualitative characteri- 
zation of the entanglement for arbitrary (pure or mixed) two- 
modes Gaussian states. Note that v_ takes a particularly sim- 
ple form for entangled symmetric states, whose standard form 
has a = b 



(38) 



Inserting Eqs. Jl 5i into Eq. d37l yields the following condi- 
tion for a two-mode squeezed thermal state £ Vi r to be entan- 
gled 



smli (2r) > 



(39) 



As for the quantification of entanglement, no fully satisfac- 
tory measure is known at present for arbitrary mixed bipartite 
Gaussian states. However, a quantification of entanglement 
which can be computed for general Gaussian states is pro- 
vided by the negativity Af, first introduced in Ref. 13911 - later 
thoroughly discussed and extended in Ref. 1 24] to CV sys- 
tems. The negativity of a quantum state g is defined as 



Af(g) 



(40) 



where g is the partially transposed density matrix and ||o|| 1 = 
Tr|d| stands for the Schatten 1-norm (the so called 'trace 
norm') of the hermitian operator 6. The quantity Af(g) is 
equal to IX^-M' me modulus of the sum of the negative 
eigenvalues of g, quantifying the extent to which g fails to 
be positive. Strictly related to Af is the logarithmic negativity 
Erf, defined as Ejy = ln||g||i, which constitutes an upper 
bound to the distillable entanglement of the quantum state g. 

The negativity has been proven to be convex and monotone 
under LOCC (local operations and classical communications) 
J40I1 . but fails to be continuous in trace norm on infinite di- 
mensional Hilbert spaces. However, this problem can be cir- 
cumvented by restricting to physical states with finite mean 
energy 1 14j] . 

We will now show that for any two-mode Gaussian state 
g the negativity is a simple decreasing function of z>_, which 
is thus itself an entanglement monotone. The norm || r ||i is 
unitarily invariant; in particular, it is invariant under global 
symplectic operations in phase space. Considering the sym- 
plectic diagonalization v = diag (£_, z>_, D+, v + ) of <r, this 
means that Af(g) = (||£||i — l)/2- Now, because of the mul- 
tiplicativity of the norm || • ||i, we have just to compute the 
trace of the single-mode thermal-like operator p®. 



00 , _ 1 



\k)(k\ 



Such an operator is obviously normalized for vi > 1, yielding 
\\p®. ||i = 1 (all eigenvalues are positive). Instead, if v.i < 1, 
then || pf. || 1 = 1/ui. The proof is completed by showing that 
the largest symplectic eigenvalue v + of e> fulfills v + > 1. 
This is obviously true for separable states, therefore we can 
set c + c_ < 0. With this position, it is easy to show that 

i>+ > > 1 . 



Thus, we obtain \\pf ||i = 1 and 



0i 



J\f(g) — max 



0, 



1 - i>_ 



2/>_ 



En{q) — max [0, — lnr^ 



(41) 



(42) 



This is a decreasing function of the smallest partially trans- 
posed symplectic eigenvalue £>_, quantifying the amount by 
which inequality d36i is violated. The eigenvalue i>_ thus 
completely qualifies and quantifies the quantum entanglement 
of a Gaussian state er. Notice that, in such an instance, this 
feature holds for nonsymmetric states as well. 

We finally mention that, as far as symmetric states are con- 
cerned, another measure of entanglement, the entanglement 
of formation Ep 14 IB . can be actually computed 1 28 ] . Since 
Ep turns out to be, again, a decreasing function of it pro- 
vides for symmetric states a quantification of entanglement 
fully equivalent to the one provided by the logarithmic nega- 
tivity Epj. 
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SYMPLECTIC EIGENVALUES AND EPR 
CORRELATIONS 



such a constraint, minimizing £ Vi becomes a trivial matter and 
yields 



A deeper insight on the relationship between correlations 
and the eigenvalue £>_ is provided by the following observa- 
tion, which holds in the symmetric case. 

Let us define the EPR correlation £ 1 28, 29] of a continuous 
variable two-mode quantum state as 



Trer 



013 + C24 



(43) 



where 65 



for an operator 6. If £ > 1 then the 



state does not possess nonlocal correlations . The idealized 
EPR-like state 1 30] (simultaneous eigenstate of the commut- 
ing observables x\ — £2 and pi + p 2 ) has £ = 0. As for 
standard form states, one has 



Sr, 



1+P2 



= a + b - 2c + , 
= a + b + 2c_ , 



£ = a + b — 



(44) 
(45) 
(46) 



Notice that £ is not by itself a good measure of correlation be- 
casue, as one can easily verify, it is not invariant under local 
symplectic operations. In particular, applying local squeez- 
ings with parameters r, = In tii and local rotations with angles 
ifii to a standard form state, we obtain 



a ' 2 , 1 n , h i 2 1 



Zvi,e = ~(vi + -o) + ^{ v 2 + -^)-{c+v 1 v 2 - 
with 1? = (pi + if2- Now, the quantity 



viv 2 



-) COS# , 

(47) 



£ = min£„ 4 * 

has to be Sp^M) © Sp(2M) invariant. It corresponds to 
the maximal amount of EPR correlations which can be dis- 
tributed in a two-mode Gaussian state by means of local op- 
erations. Minimization in terms of § is immediate, yielding 

£ = min Wj £„ i , with 



a 1 2 1, 
2^ + ^ 



& 2 1 C_ 

-{v 2 + —)- \C + ViV 2 

2 vi viv 2 



The gradient of such a quantity is null if and only if 



1 



— - \c+\viv 2 - 



v\ j ] r |r,/' 2 

Vn 



V1V2 

viv 2 



= 0, 
= 0, 



■ (48) 

(49) 
(50) 



where we introduced the position c + c_ < 0, necessary to 
have entanglement. Eqs. ( 1491 15 0> can be combined to get 



1 



1 



(51) 



£ = 2V(a-|c+|)(a-|c_|) = 2*>_ 



(52) 



We thus see that the smallest symplectic eigenvalue of the par- 
tially transposed state is endowed with a direct physical inter- 
pretation: it quantifies the greatest amount of EPR correla- 
tions which can be created in a Gaussian state by means of 
local operations. 

As can be easily shown by a numerical analysis, such a sim- 
ple interpretation is lost for nonsymmetric states. This fact 
properly exemplifies the difficulties of handling optimization 
problems in nonsymmetric instances, encountered, e.g. in the 
computation of the entanglement of formation of such states. 



VI. PARAMETRIZATION OF GAUSSIAN STATES WITH 
SYMPLECTIC INVARIANTS 

Two-mode Gaussian states can be classified according to 
the values of their four Sp^u) ® <Sp(2,R) invariants a, 6, 
c+ and c_, which determine their standard form. It is rel- 
evant to provide a reparametrization of standard form states 
in terms of invariants which admit a direct interpretation for 
generic Gaussian states. Such invariants will be the global pu- 
rity fi = Tr q 2 , the marginal purities of the reduced states 
Hi = Trj[( Tij^ig) 2 ] and the Spu^) invariant A, whose 
meaning will become soon clear. For a two-mode Gaussian 
state one has 



/'1 



A 



(53) 



1 1 

~ . M2 = T , 

a b 
Deter = {ab) 2 - ab(c\ + (?_) + (c+c_) 2 , (54) 
2 , b 2 



a 



2c + c 



(55) 



Eqs. J53I55I > can be inverted to provide the following 
parametrization 



1 

Mi 



1 

M2 



(56) 



\ 



M1M2 



A 



(Mi - 
MIM2 



± e 



(57) 



with e = — 
4 



1 /[(Ml + M2) 2 -M1M2A] 4/zi/i2 



M1M2 



The global and marginal purities range from to 1, con- 
strained by the condition 



M > M1M2 



(58) 



Restricting to the symmetric (a = b) entangled (=$■ c+C- < 0) 
case, Eq. (15 1 1 and the fact that Vi > imply vi = v%. Under 



that can be easily shown to be a direct consequence of inequal- 
ity (^J. Notice that inequality J58I entails that no Gaussian 
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LPTP (less pure than product) states exist, at variance with the 
case of two-qubit systems |6]. 

The smallest symplectic eigenvalues of the covariance ma- 
trix a and of its partial transpose <r are promptly determined 
in terms of symplectic invariants 



2vi 



2v 



r/2 



A 



(59) 



where A = -A + 2//x| + 

The parametrization provided by Eqs. d56l I57> describes 
physical states if the radicals in Eqs. j57l I59> exist and the 
Heisenberg uncertainty principle is satisfied. All these 
conditions can be combined and recast as upper and lower 
bounds on the global symplectic invariant A 



2 + (pi-^2) 2 
/i 



2 2 



< min 



(Mi 



M2J 



< A 



2 1 
--, 1 + — 



(60) 



Let us investigate the role played by the invariant A in the 
characterization of the properties of Gaussian states. To this 
aim, we just analyse the dependence of the eigenvalue on 
A 



dv 



r,2 



d A 




1 > 0. 



(61) 



The smallest symplectic eigenvalue of the partially transposed 
state is strictly monotone in A. Therefore the entanglement of 
a generic Gaussian state cr with global purity /i and marginal 
purities /ii^ strictly increases with decreasing A. The invari- 
ant A is thus endowed with a direct physical interpretation: at 
given global and marginal purities, it determines the amount 
of entanglement of the state. 

Because, due to inequality ( I60> . A possess both lower and 
upper bounds, not only maximally but also minimally entan- 
gled Gaussian states exist. This elucidates the relations be- 
tween the entanglement and the purity of two-mode Gaussian 
states: the entanglement of such states is tightly bound by the 
amount of global and marginal purities, with only one remain- 
ing degree of freedom related to the invariant A. 



VII. EXTREMAL ENTANGLEMENT AT FIXED GLOBAL 
AND LOCAL PURITIES 

We now aim to characterize extremal (maximally and mini- 
mally) entangled Gaussian states for fixed global and marginal 
purities. As it is clear from Eq. J53l >. the standard form of 
states with fixed marginal purities always satisfies a = 1/fii, 
b = 1//U2- Therefore the complete characterization of maxi- 
mally and minimally entangled states is achieved by specify- 
ing the expression of their standard form coefficients c T . 

Let us first consider the states saturating the lower bound 
in Eq. J60i . which entails maximal entanglement. They are 



Gaussian maximally entangled states for fixed global and lo- 
cal purities (GMEMS), admitting the following parametriza- 
tion 



c± = ± 



1 



1 



(62) 



It is easily seen that such states belong to the class of asym- 
metric two-mode squeezed thermal states, with squeezing pa- 
rameter and symplectic spectrum given by 



tanh2r = 2(^i/i 2 - nlnl/ 1 ''/ '{^ + M2) , 



(63) 



1 

/' 



(Ml - _ |A*1 - M2I /Ol-i«2) 2 



2 2 

MIM2 



(64) 

In particular, any GMEMS can be written as an entangled two- 
mode squeezed thermal states [satisfying Ineq. J39il. This 
provides a characterization of two-mode thermal squeezed 
states as maximally entangled states for given global and 
marginal purities. We can restate this result as follows: given 
an initial tensor product of (generally different) thermal states, 
the unitary operation providing the maximal entanglement for 
given values of the local purities /i/s is given by a two-mode 
squeezing, with squeezing parameter determined by Eq. 1631 . 
Notice that fixing the values of the local purities is necessary 
to attain maximal, finite entanglment; in fact, fixing only the 
value of the global purity [i always allows for an arbitrary 
large logarithmic negativity (as v_ goes to 0), achievable by 
means of two-mode squeezing with arbitrary large squeezing 
parameter. 

Nonsymmetric two-mode thermal squeezed states turn out 
to be separable in the range 



< 



A*lM2 



/il + // 2 - MlM2 



(65) 



In such a separable region in the space of purities, no entan- 
glement can occur for states of the form of Eq. 1621 . while, 
outside this region, they are properly GMEMS. We remark 
that, as a consequence, all Gaussian states whose purities fall 
in the separable region defined by inequality j65t are separa- 
ble. 

We now consider the states that saturate the upper bound in 
Eq. ( I60> . They determine the class of Gaussian least entangled 
states for given global and local purities (GLEMS). Violation 
of inequality d65l implies that 



1 

1 + — 

M" 



< 



2 2 

WM2 



2 



Therefore, outside the separable region, GLEMS fulfill 

1 



A = l + - 



(66) 



Eq. i66\ expresses partial saturation of Heisenberg relation 
Eq. H2i . Namely, considering the symplectic diagonalization 
of Gaussian states and Eq. H41 . it immediately follows that 
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FIG. 3: Plot of the maximal logarithmic negativity of Gaussian states 
for fixed marginal purities. The surface represents GMEMMS, states 
that saturate inequality J70L 



the Spurn invariant condition d66i is fulfilled if and only if 
the symplectic spectrum of the state takes the form i/_ = 1, 
v+ = 1/fi. We thus find that GLEMS are characterized by a 
peculiar spectrum, with all the mixedness concentrated in one 
'decoupled' quadrature. Moreover, by Eqs. ( I66> and dl 2i it 
follows that GLEMS are minimum-uncertainty mixed Gaus- 
sian states. They are determined by the standard form corre- 
lation coefficients 




Quite remarkably, following the analysis presented in Sec.lIIII 
it turns out that the GLEMS at fixed global and marginal pu- 
rities are also states of minimal global p— entropy for p < 2, 
and of maximal global p— entropy for p > 2. 

According to the PPT criterion, GLEMS are separable only 



if fi < /x^a/VMi + M2 

MlM2 



M?mI' Therefore, in the range 



< p < 



M1M2 



14 



? 9 

MIM2 



(68) 



Ml + M2 - MlM2 y/J4 

both separable and entangled states can be found. Instead, the 
region 

M1M2 



p > 



VJ4 



2 

M2 



9 2 

MIM5 



(69) 



can only accomodate entangled states. 

The very narrow region defined by inequality d68l > is thus 
the only region of coexistence of both entangled and sepa- 
rable Gaussian mixed states. We mention that the sufficient 
condition for entanglement ( I69l l. first derived in Ref. 1 13] has 
been independently rediscovered in another recent work 11711 . 

The upper and lower inequalities J60l > lead to the following 
further constraint between the global and the local purities 



M< 



M1M2 



M1M2 + I Mi - M2| 



(70) 



Degrees of Purity 


Separability 


p. < P1P2 


unphysical region 


P1P2 < p < — piiiia 

— — Ml +M2 _ Ml M2 


separable states 


M1M2 ^ ., ^ M1M2 


coexistence region 


M!+M2-M1M2 " h ~ ^2 +fl 2 _ ^2 p 2 


M1M2 ^ IX ^ M1M2 


entangled states 


V^f+Ml-M?^ " ~ M1M2 + IM1-M2 


U > M1 , M2 

^ M1M2 + IM1-M2 


unphysical region 



TABLE I: Classification of two-mode Gaussian states and of their 
properties of separability according to their degrees of global purity 
fi and of marginal purities ^ and ji2- 



For purities which saturate inequality (17 0> . one obtains the 
states of maximal global purity with given marginal purities. 
In such an instance GMEMS and GLEMS coincide and we 
have a unique class of entangled states depending only on the 
marginal purities Mi, 2- They are Gaussian maximally entan- 
gled states for fixed marginals (GMEMMS). In Fig.|3]the log- 
arithmic negativity of GMEMMS is plotted as a function of 
Mi,2> showing how the maximal entanglement achievable by 
Gaussian states rapidly decreases with increasing difference 
of marginal purities, in analogy with finite-dimensional sys- 
tems For symmetric states (/ii — /i 2 ) inequality ( 17 0> 
reduces to the trivial bound fi < 1 and GMEMMS reduce to 
pure two-mode squeezed vacua. 

The previous necessary or sufficient conditions for entan- 
glement are collected in Table[l]and allow a graphical display 
of the behavior of the entanglement of mixed Gaussian states 
as shown in Fig.|4] These relations classify the properties of 
separability of all two-mode Gaussian states according to their 
degree of global and marginal purities. 



A. Realization of extremally entangled states in experimental 
settings 

As we have seen, GMEMS are two-mode squeezed ther- 
mal states, whose general covariance matrix is described by 
Eqs. (|5) and d!5l >. A realistic instance giving rise to such states 
is provided by the dissipative evolution of an initially pure 
two-mode squeezed vacuum created, e.g., in a non degener- 
ate parametric down conversion process. Let us denote by <r r 
the covariance matrix of a two mode squeezed vacuum with 
squeezing parameter r, derived from Eqs. dl5> with v T = 1. 
The interaction of this initial state with a thermal noise results 
in the following dynamical map describing the time evolution 
of the covariance matrix cr(t) I3lll 



-rt. 



(1 



-rts 



(71) 



where T is the coupling to the noisy reservoir (equal to the 



inverse of the damping time) and cr 7 



i=V 



ito is the 



covariance matrix of the thermal noise. The average number 
of thermal photons rti given by 



1 



n, = 



exp (huji/kBT) — 1 



10 




FIG. 4: Summary of entanglement properties of two-mode (nonsymmetric) Gaussian states in the space of marginal purities /ii^ (x- and 
y-axes) and global purity /j,. In fact, on the z-axis we plot the ratio /i//ii/i2 to gain a better graphical distinction between the various regions. 
In this space, all physical states lay between the horizontal plane z = 1 representing product states, and the upper limiting surface representing 
GMEMMS. Separable and entangled states are well separated except for a narrow region of coexistence (depicted in black). Separable states 
fill the region depicted in dark grey, while in the region containing only entangled states we have depicted the average logarithmic negativity 
Eq. <88> . growing from white to medium grey. The mathematical relations defining the boundaries between all these regions are collected in 
TableU The three-dimensional envelope is cut at z = 3.5. 



in terms of the frequencies of the modes and of the temper- 
ature of the reservoir T. It can be easily verified that the co- 
variance matrix Eq. ( 17 11 defines a two-mode thermal squeezed 
state, generally nonsymmetric (for m ^ 712)- However, no- 
tice that the entanglement of such a state cannot persist in- 
definetely, becaus after a given time inequality d39l will be 
violated and the state will evolve into a non entangled two- 
mode squeezed thermal state. We also notice that the relevant 
instance of pure loss {n\ = rii = 0) allows the realization of 
symmetric GMEMS. 

Concerning the experimental characterization of minimally 
entangled states, one can envisage several explicit experimen- 
tal settings for their realiazion. For instance, let us consider 
(see Fig. |3 a beam splitter with transmittivity r = 1/2 (corre- 
sponding to a two-mode rotation of angle 7r/4 in phase space). 
Suppose that a single-mode squeezed state, with covariance 



that reads 



diag ( 



■2r 



-2r\ 



(like, e.g., the result of a de- 



matnx <T lr 

generate parametric down conversion in a nonlinear crystal), 
enters in the first input of the beam splitter. Let the other in- 
put be an incoherent thermal state produced from a source at 
equilibrium at a temperature T. The purity \i of such a state 
can be easily computed in terms of the temperature T and of 
the frequency of the thermal mode u> 

exp (Sw/igT) — 1 



1 



(72) 



exp (hui/ksT) 

The state at the output of the beam splitter will be a corre- 
lated two-mode Gaussian state with covariance matrix <r out 



I n + k 







n — k 



& out 2 
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OPO / OPA 
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FIG. 5: Possible scheme for the generation of Gaussian least en- 
tangled mixed states (GLEMS). A single-mode squeezed state (ob- 
tained, for example, by an optical parametric oscillator or amplifier) 
interferes with a thermal state through a 50:50 beam splitter. The 
resulting two-mode state is a minimally entangled mixed Gaussian 
state at given global and marginal purities. 
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with k = e 2r and n — /i~ 1 . By immediate inspection, the 
symplectic spectrum of this covariance matrix is v_ = 1 and 
v+ = Therefore the output state is always a state with 
extremal generalized entropy at a given purity (the state can 
be seen as the tensor product of a vacuum and of a thermal 
state, on which one has applied symplectic transformation). 
Moreover, the state is entangled if 



cosh(2r) > 



1 _ exp(2Huj/k B T) + 1 



2fi exp(2Sw/fc B T) - 1 



(73) 



Tuning the experimental parameters to meet the above con- 
dition, indeed makes the output state of the beam splitter a 
symmetric GLEMS. It is interesting to observe that nonsysm- 
metric GLEMS can be produced as well by choosing a beam 
splitter with transmissivity different from 0.5. 



VIII. EXTREMAL ENTANGLEMENT AT FIXED GLOBAL 
AND LOCAL GENERALIZED ENTROPIES 



In this section we introduce a more general characterization 
of the entanglement of generic two-mode Gaussian states, by 
exploiting the generalized p— entropies as measures of global 
and marginal mixedness. For ease of comparison we will 
carry out this analysis along the same lines followed in the 
previous Section, by studying the explicit behavior of the 
global invariant A, directly related to the logarithmic nega- 
tivity Ej^f at fixed global and marginal entropies. This study 
will clarify the relation between A and the generalized en- 
tropies S p and the ensuing consequences for the entanglement 
of Gaussian states. 

We begin by observing that the standard form covari- 
ance matrix <x of a generic two-mode Gaussian state can be 
parametrized by the following quantities: the two marginals 
/j,i } 2 (or any other marginal S Pl 2 because all the local, single- 
mode entropies are equivalent for any value of the integer p), 
the global p— entropy S p (for some chosen value of the integer 
p), and the global symplectic invariant A. On the other hand, 
Eqs. (I18I21I13> provide an explicit expression for any S p as a 
function of p and A. Such an expression can be exploited to 
study the behaviour of A as a function of the global purity /i, 
at fixed marginals and global S p (from now on we will omit 
the explicit reference to fixed marginals). One has 



OA 



2_ OR 
R 2 dA 



2 dS p /dA\ 



R 



Sp R 2 dS p /dR\ A 
2 N P (A,R) 



R 2 D P (A,R) ' 
where we have defined the inverse participation ratio 

2 



(74) 



R 



(75) 



and the remaining quantities N p and D p read 

N P (A,R) 



(R + 2 + 2VA + R) p - 1 
(R + 2 - 2VATi?) p - 1 
(R - 2 + 2y/A~R) p - 1 



VA~R 
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- R- 
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2 + 2VA - 


-Rf- 1 



y/A + R 



(76) 



Now, it is easily shown that the ratio N p (A, R) / D p (A, R) is 
increasing with increasing p and has a zero at p = 2 for any 
A, R; in particular, its absolute minimum (—1) is reached in 
the limit (A — * 2, R — ► 2, p — > 1). Thus the derivative 
Eq. (1741 is negative for p < 2, null for p = 2 (in this case A 
and S*2 = 1 — /i are of course regarded as independent vari- 
ables) and positive for p > 2. This implies that, for given 
marginals, keeping fixed any global S p for p < 2 the min- 
imum (maximum) value of A corresponds to the maximum 
(minimum) value of the global purity /i. Instead, by keeping 
fixed any global S p for p > 2 the minimum of A is always 
attained at the minimum of the global purity /i. 

This observation allows to determine rather straightfor- 
wardly the states with extremal A. They are extremally entan- 
gled states because, for fixed global and marginal entropies, 
the logarithmic negativity of a state is determined only by the 
one remaining independent global symplectic invariant, rep- 
resented by A in our choice of parametrization. If, for the 
moment being, we neglect the fixed local purities, then the 
states with maximal A are the states with minimal (maximal) 
[i for a given global S p with p < 2 (p > 2) (see Sec. |lll| 
and Fig. [Q. As found in Sec. [HI] such states are minimum 
uncertainty two-mode states with mixedness concentrated in 
one quadrature. We have shown in Sec. IV11I that they cor- 
respond to Gaussian least entangled mixed states (GLEMS) 
whose standard form is given by Eq. (I67> . As can be seen 
from Eq. ( I67L these states are consistent with any legitimate 
physical value of the local invariants /ii,2- We therefore con- 
clude that all Gaussian states with maximal A for any fixed 
triple of values of global and marginal entropies are GLEMS. 

Viceversa one can show that all Gaussian states with min- 
imal A for any fixed triple of values of global and marginal 
entropies are Gaussian maximally entangled mixed states 
(GMEMS). This fact is immediately evident in the symmetric 
case because the extremal surface in the S p vs. Sl diagrams is 
always represented by symmetric two-mode squeezed states 
(symmetric GMEMS). These states are characterized by a 
degenerate symplectic spectrum and encompass only equal 
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choices of the local invariants: \i\ = p,2- In the nonsymmetric 
case, the given values of the local entropies are different, and 
the extremal value of A is further constrained by inequality 

(EDI 



A-R> 



(Mi ~ M2) s 

? 2 

MIM2 



From Eq. ( 1741 it follows that 



d(A - R) 



dA 



N P (A,R) 
D P (A,R) 



> 0. 



(77) 



(78) 



because N p (A, R)/D p (A, R) > -1. Thus, A - R is an in- 
creasing function of A at fixed /X1.2 and S p , and the minimal 
A corresponds to the minimum of A — R, which occurs if 
inequality 1771 is saturated. Therefore, also in the nonsym- 
metric case, the two-mode Gaussian states with minimal A at 
fixed global and marginal entropies are GMEMS. 

Summing up, we have shown that GMEMS and GLEMS in- 
troduced in the previous section at fixed global and marginal 
linear entropies, are always extremally entangled Gaussian 
states, whatever triple of generalized global and marginal en- 
tropic measures one chooses to fix. Maximally and minimally 
entangled states of CV systems are thus very robust with re- 
spect to the choice of different measures of mixedness. This 
is at striking variance with the case of discrete variable sys- 
tems, where it has been shown that fixing different measures 
of mixedness yields different classes of maximally entangled 
states 01 . 

Furthermore, we will now show that the characterization 
provided by the generalized entropies leads to some remark- 
able new insight on the behavior of the entanglement of CV 
systems. The crucial observation is that for a generic p, the 
smallest symplectic eigenvalue of the partially transposed co- 
variance matrix, at fixed global and marginal entropies, is 
not in general a monotone function of A, so that the connec- 
tion between extremal A and extremal entanglement turns out 
to be, in some cases, inverted. In particular, while for p < 2 
the GMEMS and GLEMS surfaces tend to be more separated 
as p decreases, for p > 2 the two classes of extremally en- 
tangled states get closer with increasing p and, within a par- 
ticular range of global and marginal entropies, they exchange 
their role. GMEMS (i.e. states with minimal A) become min- 
imally entangled states and GLEMS (i.e. states with maximal 
A) become maximally entangled states. This inversion always 
occurs for all p > 2. 

To understand this interesting behaviour, let us study the 
dependence of the symplectic eigenvalue V- on the global in- 
variant A at fixed marginals and at fixed S p for a generic p. 
Using Maxwell's relations, we can write 



d(29 2 _) 



dA 



Clearly, for 



d(2D 2 _) 



d(2vl) 8S p /dA\ R 
dA R ' OR A ' dS p /dR\ A ' 

(79) 

> GMEMS and GLEMS retain their usual 
interpretation, whereas for n p < they exchange their role. 



On the node k v 



GMEMS and GLEMS share the same 
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FIG. 6: Plot of the nodal surface which solves the equation k p = 
with k p defined by Eq. 1801 . for (a) p = 3 and (b) p = 4. The 
entanglement of Gaussian states that lie on the leaf-shaped surfaces 
is fully quantified in terms of the marginal purities and the global 
generalized entropy (a) S3 or (b) S4. The equations of the surfaces 
in the space £ p = {/j,i,/j,2, S p } are given by Eqs. !83H85i . 



at k p — is fully determined by the global and marginal 
p— entropies alone, and does not depend any more on A. Such 
nodes also exist in the case p < 2 in two limiting instances: 
in the special case of GMEMMS (states with maximal global 
purity at fixed marginals) and in the limit of zero marginal 
purities. We will now show that, besides these two limiting 
behaviors, a nontrivial node appears for all p > 2, implying 
that on the two sides of the node GMEMS and GLEMS in- 
deed exhibit opposite behaviors. Because of Eq. i74l , k p can 
be written in the following form 



Kp = K 2 



R 



N P (A,R) 
VA 2 - R 2 D P (A,R) 



(80) 



with N p and D p defined by Eq. d76l > and 



K 2 = -1 + 



entanglement, i.e. the entanglement of all Gaussian states 



9 9 ' 
Mi M2 

A 

VA 2 - R 2 ' 



The quantity k p in Eq. J80l > is a function of p, R, A, and 
of the marginals; since we are looking for the node (where 
the entanglement is independent of A), we can investigate the 
existence of a nontrivial solution to the equation k p = fixing 
any value of A. Let us choose A = 1 +i? 2 /4 that saturates the 
Heisenberg uncertainty relation and is satisfied by GLEMS. 
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With this position, Eq. J80l > becomes 

R 



(81) 



The existence of the node depends then on the behavior of the 
function 



fp(R) 



2 [(R + 2)p~ 2 - (R - 2)p~ 2 ] 
(R + 4)(i? + 2)P" 2 - (R - 4)(i? - 2) 



(82) 



In fact, as we have already pointed out, k 2 is always posi- 
tive, while the function f p (R) is an increasing function of p 
and, in particular, it is negative for p < 2, null for p = 2 
and positive for p > 2, reaching its asymptote 2/(R + 4) in 
the limit p — > oo. This entails that, for p < 2, k p is always 
positive, which in turn implies that GMEMS and GLEMS are 
respectively maximally and minimally entangled two-mode 
states with fixed marginal and global p— entropies in the range 
p < 2 (including both von Neumann and linear entropies). 
On the other hand, for any p > 2 one node can be found 
solving the equation ft p (/ii, fi2, 2/m) = 0. The solutions 
to this equation can be found analytically for low p and nu- 
merically for any p. They form a continuum in the space 
{fii, fi2, m} which can be expressed as a surface of general 
equation /i = ^(/zi,/i2)- Since the fixed variable is S p and 
not fi it is convenient to rewrite the equation of this surface 
in the space £ p = {/zi, /j,2, S p }, keeping in mind the rela- 
tion ( I34t . holding for GLEMS, between [i and S p . In this way 
the nodal surface (k p — 0) can be written in the form 



S p = Sp(pi,fi2) 



1 



9 P [(Mp(A*i>/"2)) x ] 



(83) 



The entanglement of all Gaussian states whose entropies lie 
on the surface S£(px, (12) is completely determinated by the 
knowledge of fix, /i 2 and S p . The explicit expression of the 
function n p (fii, ^2) depends on p but, being the global purity 
of physical states, is constrained by the inequality 



M1M2 < /ip(Mi,M2) < 



MlM2 



/il/i2 + |Mi - 

The nodal surface of Eq. J83l > constitutes a 'leaf, with base at 
the point /z£ (0,0) = and tip at the point //£(a/3/2, a/3/2) 
1 . for any p > 2; such a leaf becomes larger and flatter with 
increasing p (see Fig.|6j. Forp > 2, the function f p (R) de- 
fined by Eq. d82l is negative but decreasing with increasing R, 
that is with decreasing /i. This means that, in the space of en- 
tropies S p , above the leaf (S p > Sg) GMEMS (GLEMS) are 
still maximally (minimally) entangled states for fixed global 
and marginal generalized entropies, while below the leaf they 
are inverted. Notice also that for /11.2 > Vo/2 no node and 
so no inversion can occur for any p. Each point on the leaf- 
shaped surface of Eq. J83I corresponds to an entire class of in- 
finitely many two-mode Gaussian states (including GMEMS 



and GLEMS) with the same marginals and the same global 
S p = Sp (/ii, fj,2), which are all equally entangled, since their 
logarithmic negativity is completely determined by /ii, //2 and 
S p . For the sake of clarity we provide the explicit expressions 
of Mp(Mi) A*2)> as plotted in Fig.|6]for the cases (a) p = 3, and 

(b)p = 4. 



^1,^2) = 3 3 , (84) 

M4(l»l>M2) = V3M1M2 / ((4 + A ~ 2 MlM2 + 

\J (m? + 14) (m? + i4- (4i4) + 1A1A 

(85) 

Apart from the relevant 'inversion' feature shown by 
p— entropies for p > 2, the possibility of an accurate charac- 
terization of CV entanglement based on global and marginal 
entropic measures still holds in the general case for any p. In 
particular, the set of all Gaussian states can be again divided, 
in the space of global and marginal S p 's, into three main ar- 
eas: separable, entangled and coexistence region. It can be 
thus very interesting to investigate how the different entropic 
measures chosen to quantify the degree of global mixedness 
(all marginal measures are equivalent) behave in classifying 
the separability properties of Gaussian states. Fig.0provides 
a numerical comparison of the different characterizations of 
entanglement obtained by the use of different p— entropies, 
with p ranging from 1 to 4, for symmetric Gaussian states 
(S Pl = S P2 = S Pi ). The last restriction has been imposed just 
for ease of graphical display. The following considerations, 
based on the exact numerical solutions of the transcendental 
conditions, will take into account nonsymmetric states as well. 

The mathematical relations expressing the boundaries be- 
tween the different regions in Fig. [7] are easily obtained for 
any p by starting from the relations holding for p = 2 (see Ta- 
ble Q and by evaluating the corresponding for each 
1^(1^1,2)- For any physical symmetric state such a calculation 
yields 



0< (p-l)S p < l-g P 

=> entangled, 



Mi 



<(p-l)Sp < l-g 2 p ' 2 >' 



coexistence, 



separable. (86) 

Equations d86l > were obtained exploiting the multiplicativity 
of p— norms on product states and using Eq. (I27t for the lower 
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FIG. 7: Summary of the entanglement properties for symmetric Gaussian states at fixed global and marginal generalized p— entropies, for (a) 
p = 1 (von Neumann entropies), (b) p = 2 (linear entropies), (c) p = 3, and (d) p = 4. In the entangled region, the average logarithmic 
negativity Ej^(S Pi , S p ) Eq. J88I is depicted, growing from gray to black. For p > 2 an additional dashed curve is plotted; it represents the 
nodal line of inversion. Along it the entanglement is fully determined by the knowledge of the global and marginal generalized entropies 
S Pi ,S p , and GMEMS and GLEMS are equally entangled. On the left side of the nodal line GMEMS (GLEMS) are maximally (minimally) 
entangled Gaussian states at fixed S Pi , S p . On the right side of the nodal line they are inverted: GMEMS (GLEMS) are minimally (maximally) 
entangled states. Also notice how the dashed region of coexistence becomes narrower with increasing p. The equations of all boundary curves 
can be found in Eq. 1861 . 



boundary of the coexistence region (which represents GLEMS 
becoming entangled) and Eq. d30i for the upper one (which 
expresses GMEMS becoming separable). Let us mention also 
that the relation between any local entropic measure S Pi and 
the local purity //$ is obtained directly from Eq. i !21t and reads 



= 1 - gp(lM) (87) 
p-1 

We notice prima facie that, with increasing p, the entan- 
glement is more sharply qualified in terms of the global and 
marginal p— entropies. In fact the region of coexistence be- 
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tween separable and entangled states becomes narrower with 
higher p. Thus, somehow paradoxically, with increasing p the 
entropy S p provides less information about a quantum state, 
but at the same time it yields a more accurate characterization 
and quantification of its entanglement. In the limit p — > oo 
all the physical states collapse to one point at the origin of the 
axes in the space of generalized entropies, due to the fact that 
the measure Sao is identically zero. 



IX. QUANTIFYING ENTANGLEMENT: THE AVERAGE 
LOGARITHMIC NEGATIVITY 



We have extensively shown that knowledge of the global 
and marginal generalized p— entropies accurately character- 
izes the entanglement of Gaussian states, providing strong 
sufficient or necessary conditions. The present analysis nat- 
urally leads us to propose an actual quantification of entan- 
glement, based exclusively on marginal and global entropic 
measures, enriching and generalizing the approach introduced 
inRef. CI- 

Outside the separable region, we can formally define the 
maximal entanglement Ejyf max (S Pl 2 , S p ) as the logarithmic 
negativity attained by GMEMS (or GLEMS, below the inver- 
sion nodal surface for p > 2, see Fig.|S}. In a similar way, in 
the entangled region GLEMS (or GMEMS, below the inver- 
sion nodal surface for p > 2) achieve the minimal logarithmic 
negativity Ej^f max {S Pl 2 , S p ). The explicit analytical expres- 
sions of these quantities are unavailable for any p / 2 due 
to the transcendence of the conditions relating S p to the sym- 
plectic eigenvalues. The surfaces of maximal and minimal en- 
tanglement in the space of the global and local S p are plotted 
in Fig. [8] for symmetric states. In the plane S p = the up- 
per and lower bounds correctly coincide, since for pure states 
the entanglement is completely quantified by the marginal en- 
tropy. For mixed states this is not the case but, as the plot 
shows, knowledge of the global and marginal entropies strictly 
bounds the entanglement both from above and from below. 
For p > 2, we notice how GMEMS and GLEMS exchange 
their role beyond a specific curve in the space of S p 's. The 
equation of this nodal curve is obtained from the general leaf- 
shaped nodal surfaces of Eqs. J83H85I . by imposing the sym- 
metry constraint (S Pl = S P2 = S Pi ). We notice again how the 
S p 's with higher p provide a better characterization of the en- 
tanglement, even quantitatively. In fact, the gap between the 
two extremally entangled surfaces in the S p 's space becomes 
smaller with higher p. Of course the gap is exactly zero all 
along the nodal line of inversion for p > 2. 

We will now introduce a particularly convenient quantita- 
tive estimate of the entanglement based only on the knowl- 
edge of the global and marginal entropies. Let us define the 
"average logarithmic negativity" E_\r as 



Gaussian states. To this aim, we define the relative error <5EV 
on Ejj as 



Eu(S Pl 2 , S p ) 



Em max(Sp 12 , Sp) -\- Ej\f rn i n (S pi 2 , Sp) 



6Etf(S Pl:2 , S p ) 



E N max i,Sp 1 2 , Sp) Ej^jnin {Sp 12 , Sp) 



(88) 

We will now show that this quantity, fully determined by the 
global and marginal entropies, provides a reliable quantifica- 
tion of entanglement (logarithmic negativity) for two-mode 



E]\[rnax(Sp 1 2 , Sp) ~\- -E^/Vmm 2 , Sp) 

(89) 

As Fig.|9]shows, this error decreases exponentially both with 
decreasing global entropy and increasing marginal entropies, 
that is with increasing entanglement. In general the relative 
error SEj^f is 'small' for sufficiently entangled states; we will 
present more precise numerical considerations in the subcase 
p = 2. Notice that the decaying rate of the relative error is 
faster with increasing p: the average logarithmic negativity 
turns out to be a better estimate of entanglement with increas- 
ing p. For p > 2, <5£V is exactly zero on the inversion node, 
then it becomes finite again and, after reaching a local maxi- 
mum, it goes asymptotically to zero (see the insets of Fig. [5}. 

All the above considerations, obtained by an exact numeri- 
cal analysis, show that the average logarithmic negativity E^ 
at fixed global and marginal p— entropies is a very good es- 
timate of entanglement in CV systems, whose reliability im- 
proves with increasing entanglement and, surprisingly, with 
increasing order p of the entropic measures. 



A. Direct estimate of the entanglement 

In the present general framework, a peculiar role is played 
by the case p = 2, i.e. by the linear entropy Sl (or, equiv- 
alently, the purity /x). The previous general analysis on the 
whole range of generalized entropies S p , has remarkably 
stressed the privileged theoretical role of the instance p = 2, 
which discriminates between the region in which extremally 
entangled states are unambiguously characterized and the re- 
gion in which they can exchange their roles. Moreover, the 
graphical analysis shows that, in the region where no inver- 
sion takes place (p < 2), fixing the global S2 = 1 — yields 
the most stringent constraints on the logarithmic negativity of 
the states (see Figs. [8] [9J. Notice that such constraints, 
involving no transcendental functions for p = 2, can be eas- 
ily handled analytically. A crucial experimental consideration 
strengthens these theoretical and practical reasons to privilege 
the role of S2. In fact, S2 can indeed, assuming some prior 
knowledge about the state (essentially, its Gaussian character), 
be measured through conceivable direct methods, in particular 
by means of single-photon detection schemes fl7ll or of the 
upcoming quantum network architectures l42ll . No complete 
homodyne reconstruction of the density matrix is needed in 
such schemes. Very recently, a first important and promising 
step in this direction has been realized with the experimental 
implementation of direct photon detection for the measure- 
ment of the squeezing and purity of a single-mode squeezed 
vacuum with a setup that required only a tunable beam splitter 
and a single-photon detector 1 4311 . 

As already anticipated, for p — 2 we can provide analyti- 
cal expressions for the extremal entanglement in the space of 
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FIG. 8: Upper and lower bounds on the logarithmic negativity of symmetric Gaussian states as functions of the global and marginal generalized 
p— entropies, for (a) p = 1 (von Neumann entropies), (b) p = 2 (linear entropies), (c) p — 3, and (d) p = 4. The black (gray) surface represents 
GMEMS (GLEMS). Notice that forp > 2 GMEMS and GLEMS surfaces intersect along the inversion line (meaning they are equally entangled 
along that line), and beyond it they interchange their role. The equations of the inversion lines are obtained from Eqs. I83H85I . with the position 



S n 



Si) 
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global and marginal purities 

max 



- log 

2 S 



■£wWn(Ml,2)M) = -o lo S 



1 

/' 



Mi + Mg 
2/xf/il 



x \ m + (j.2 — \ (mi + M2) 2 - 



A 4 



1111 

^f + ^| _ V ~ 2 



'ii 1 1 



1 



(90) 



(91) 



Consequently, both the average logarithmic negativity SEjsf, 
defined in Eq. j88l >. and the relative error SEtf, given by 
Eq. ( I89> , can be easily evaluated in terms of the purities. The 
relative error is plotted in Fig. |9(b)| for symmetric states as a 
function of the ratio SlJSl- Notice, as already pointed out 
in the general instance of arbitrary p, how the error decays 
exponentially. In particular, it falls below 5% in the range 
Sl < (m > Mi), which excludes at most very weakly 
entangled states (states with EV < 1 1 44]). Let us remark 
that the accuracy of estimating entanglement by the average 
logarithmic negativity proves even better in the nonsymmetric 
case /ii ^ H2, essentially because the maximal allowed entan- 
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FIG. 9: The relative error 5Saa Eq. |89j on the average logarithmic negativity as a function of the ratio S Pi /S p , for (a) p = 1, (b)p = 2, (c) 
p = 3, (d) p = 4, plotted at (a) SV = 1, (b) Sl ~ 1/2, (c) S3 = 1/4, (d) S4 — 1/6. Notice how, in general, the error decays exponentially, 
and in particular faster with increasing p. For p > 2, notice how the error reaches zero on the inversion node (see the insets), then grows and 
reaches a local maximum before going back to zero asymptotically. 



glement decreases with the difference between the marginals, 
as shown in Fig. l|3} 

The above analysis proves that the average logarithmic neg- 
ativity Ej^f is a reliable estimate of the logarithmic negativ- 
ity Erf, improving as the entanglement increases. This al- 
lows for an accurate quantification of continuous variable en- 
tanglement by knowledge of the global and marginal puri- 
ties. As we already mentioned, the latter quantities may be 
in turn amenable to direct experimental determination by ex- 
ploiting recent single-photon detection proposals or the up- 
coming technology of quantum networks. 



X. SUMMARY AND CONCLUDING REMARKS 

Summarizing, we have pointed out the existence of both 
maximally and minimally entangled two-mode Gaussian 
states at fixed local and global generalized p— entropies. The 
analytical properties of such states have been studied in de- 
tail for any value of p. Remarkably, for p < 2, minimally 
entangled states are minimum uncertainty states, saturating 
the Heisenberg principle, while maximally entangled states 
are nonsymmetric two-mode squeezed thermal states. Inter- 
estingly, for p > 2 and in specific ranges of the values of 
the entropic measures, the role of such states is reversed. In 
particular, for such quantifications of the global and local en- 



tropies, two-mode squeezed thermal states, often referred to 
as continuous variable analog of maximally entangled states, 
turn out to be minimally entangled. 

In the search for the extremally entangled states, we fo- 
cused on the hierarchy of mixedness induced by p— entropies 
on the set of arbitrary multi-mode states, investigating several 
related subjects, like the ordering of such different entropic 
measures and the analytical comparison between the generic 
S p (and in particular the von Neumann entropy Sy) and the 
linear entropy Sl- Moreover, we have introduced the no- 
tion of "average logarithmic negativity" for given global and 
marginal generalized p-entropies, showing that it provides a 
reliable estimate of CV entanglement in a wide range of phys- 
ical parameters. 

Our analysis also clarifies the reasons why the linear en- 
tropy is a 'privileged' measure of mixedness in continuous 
variable systems. It is naturally normalized between and 
1, it offers an accurate qualification and quantification of en- 
tanglement of any mixed state while giving significative in- 
formation about the state itself and, crucially, is the only en- 
tropic measure which could be directly measured in the near 
future by schemes involving only single-photon detections or 
the technology of quantum networks, without requiring a full 
homodyne reconstruction of the state. 

More generally, the present analysis shows that some of the 
canonical measures of entanglement and mixedness in the dis- 
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crete variable scenario, such as the entanglement of formation 
and the von Neumann entropy, may not be the best choices 
for the characterization of mixed states of continuous variable 
systems. Discontinuous behaviors appear in the limit of in- 
finite dimensions for quantities like the entanglement of for- 
mation, even when restricted to the special class of Gaussian 
states. Entanglement measures such as the logarithmic nega- 
tivity, and informational measures such as the linear entropy 
and the generalized entropies of higher order provide quantifi- 



cations which are, in many respects, more satisfactory (both 
mathematically and physically) in the context of CV systems. 
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